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PREFACE 


Part  of  the  Project  RAND  research  program  consists  of 
basic  supporting  studies.  One  aspect  of  this  concerns  the 
solution  of  the  partial  differential  equations  of  applied 
mathematics  on  automatic  digital  computing  machines  through 
approximations  by  difference  equations. 

In  the  present  Memoreindum  the  author  investigates  the 
determination  of  discrete  approximations  that  manifestly 
possess  properties,  such  as  boundedness,  positivity,  and  so 
on,  of  the  given  differential  operator. 
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SUMMARY 

The  study  of  finite— difference  approximations  to  partial 
differential  equations  Is  well  advanced,  and  the  Increasing  use 
and  significance  of  digital  computers  guarantees  that  this  study 
will  continue  to  be  actively  pursued.  One  question,  however, 
that  does  not  seem  to  have  been  discussed  to  any  extent  Is  that 
of  finding  discrete  approximations  that  manifestly  exhibit  the 
boundedness,  positivity,  and  so  on,  of  the  original  continuous 
operator. 

Thus,  for  exaii.ple,  the  recurrence  relation 
u(t  +  a)  =  (1  —  aA)u(t),  u(0)  =  c, 
yields  the  nonnegativity  we  associate  with  the  solution  of 
u'  =  —  au,  u(0)  =  c, 

provided  that  aA  <  1.  On  the  other  hand,  the  relation 

U(x,t  +  a-g,t)  ^  u(x  - 

2 

t  =  0,A  shows  very  clearly  the  nonnegativity  associated 

with  the  solution  of 

^t  ^  ^xx'  =  g(x)* 

Approximations  of  this  type  are  very  useful  computationally 
because  of  their  stability  properties.  In  this  Memorandum,  the 
author  discusses  briefly  the  problem  of  obtaining  approximations 
of  arbitrai’ily  high  degree  with  the  same  property,  and  the 
related  problem  for 

Up  =  +  g(u). 
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SOME  QUESTIONS  CONCERNING  DIFFERENCE  APPROXIMATIONS 
TO  PARTIAL  DIFFERENTIAL  EQUATIONS 

1 .  INTRODUCTION 

The  study  of  finite— difference  approximations  to  partial 
differential  equations  Is  well  advanced,  and  the  Increasing 
use  aind  significance  of  digital  computers  guarantees  that 
this  study  will  continue  to  be  actively  pursued.  One  question, 
however,  that  does  not  seem  to  have  been  discussed  to  any 
extent  Is  that  of  finding  discrete  approximations  that 
manifestly  exhibit  the  bovuidedness,  positivity,  auid  so  on,  of 
the  original  continuous  operator. 

Thus,  for  example,  the  recurrence  relation  • 

I 

(1.1)  u(t  +  a)  «  (1  —  aA)u(t),  u(0)  ■  c, 

yields  the  nonnegativity  we  associate  with  the  solution  of 

(1.2)  u'  ■  —  au,  u(0)  «  c, 

provided  that  aA  <  1.  On  the  other  hand,  the  relation 

(1.5)  u(x.t  + 

2 

t  »  0,A  ,...,  shows  very  clearly  the  nonnegativity  associated 
with  the  solution  of 

I 

(1.4)  u^  -  u^,  u(x,0)  =.  g(x). 

Approximations  of  this  type  are  very  useful  computationally  ' 

because  of  their  stability  properties.  We  shall  discuss  I 

1 

I 
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briefly  below  the  problem  of  obtaining  approximations  of 
arbitrarily  high  degree  with  the  same  property,  amd  the 
related  problem  for 

(1.5)  =  uu^  +  g(u)* 

For  previous  work  and  computational  results,  see  [l]. 

2.  THE  ONE-DIMENSIONAL  HEAT  EQUATION 

To  obtain  higher— order  approximations  possessing  the 
property  that  nonnegativity  Is  preserved,  we  can  proceed  as 
follows.  Write 

(2.1)  u(x,t  +  A  )  =  2  w.(u(x  +  r.A,t)  +  u(x  -  r.A,t)). 

1=1  ^  ^  ^ 

Using  the  relations  u^  =  '^xx^  '^tt  “  ^xxxx' 
obtain  a  set  of  moment  relations 

(2.2)  1  =  2  w  , 

1 


and  so  on.  The  question  of  existence  of  real  r^  and 

2k\ 

positive  Wj^  yielding  am  error  which  Is  0(A  )  for  a  given 

k  and  N  can  then  be  determined  by  Invoking  classical 
moment  theory . 

More  interesting  would  be  a  general  result  to  the 
effect  that  any  linear  partial  differential  equation 
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with  constant  coefficients  can  be  approximated  by 
difference  relations  of  arbitrarily  high  degree  which  preserve 
the  nonnegativity  and  boundedness  properties  of  the  original 
equation.  Let  us  further  note  that  a  nonlinear  approximation, 
e.g.,  one  arising  from  a  branching  process,  may  be  superior  In 
tnat  It  requires  fewer  terms  to  obtain  an  equivalent 
approx Imat Ion . 

3 .  THE  EQUATION  u^  -  uu^ 

Consider  now  the  equation 

(3.1)  »  uu^,  u(x,0)  -  g(x)  , 

which  Is  often  used  (see  [1))  as  a  simple  model  of  an  equation 
generating  a  shock  wave.  Since  It  possesses  an  explicit 
analytic  solution.  It  Is  useful  In  testing  numerical  Inte— 

O 

gratlon  techniques.  An  approximation  to  terms  In  0(A  )  Is 
obtained  using  the  difference  relation 

(3.2)  u(x,t  +  A)  -  u(x  +  u(x,t)A,t), 

p 

An  approximation  to  terms  In  0(A  )  Is  obtained  from 

(3.5)  vi(x,t  +  a)  -  u(x  +  u(x  +  u(x,t)A,t)A, t) . 

One  might  suspect  that  arbitrarily  accurate  approximations 
could  be  obtained  by  continuing  In  this  fashion,  but  It  seems 
difficult  to  establish  this. 


Similarly,  the  equation 


(3.4)  =  uu^  +  g(u) 

has  the  approximation 

(3*5)  u(x,t  +  a)  =  u(x  +  u(x,t)A,t)  +  g(u)A, 
valid  to  O(A^),  and 

(3.6)  u(x,t  +  a)  =  u(x  +  u(x  +  u(x,t)A,t)A  + 

+  g(u(x  +  u(x,t)A)A  +  ^ 

valid  to  O(A^). 
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